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Method and Apparatus for Determining 
the Nature of Subterranean Reservoirs 

The present invention relates to a method and apparatus for determining the nature of 
5 submarine and subterranean reservoirs. The invention is particularly suitable for 
determining whether a reservoir, whose approximate geometry and location are 
known, contains hydrocarbons or water, though it can also be applied to detecting 
reservoirs with particular characteristics. 

10 Currently, the most widely used techniques for geological surveying, particularly in 
submarine situations, are seismic methods. These seismic techniques are capable of 
revealing the structure of the subterranean strata with some accuracy. However, 
whereas a seismic survey can reveal the location and shape of a potential reservoir, it 
cannot reveal the nature of the reservoir. 

15 

The solution therefore is to drill a borehole into the reservoir. However, the costs 
involved in drilling an exploration well tend to be in the region of £25m and since 
the success rate is generally about 1 in 10, this tends to be a very costly exercise. 

20 It is therefore an object of the invention to provide a system for determining, with 
greater certainty, the nature of a subterranean reservoir without the need to sink a 
borehole. 

r 

f 

It has been appreciated by the present applicants that while the seismic properties of 
25 oil-filled strata and water-filled strata do not differ significantly, their 
electromagnetic resistivities (permittivities) do differ. Thus, by using an 
electromagnetic surveying method, these differences can be exploited and the success 
rate in predicting the nature of a reservoir can be increased significantly. This 
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represents potentially an enormous cost saving. 

Consequently, a method and apparatus embodying these principles from the basis of 
the present applicant's co-pending International patent application PCT/GB01/00419. 

5 

This contemplates a method of determining the nature of a subterranean reservoir 
whose approximate geometry and location are known, which comprises: applying a 
time varying electromagnetic field to the strata containing the reservoir; detecting the 
electromagnetic wave field response; seeking in the wave field response, a 
component representing a refracted wave from the hydrocarbon layer; and 
determining the content of the reservoir, based on the presence or absence of a wave 
component refracted by the hydrocarbon layer. 

It also contemplates a method for searching for a hydrocarbon containing 
subterranean reservoir which comprises: applying a time varying electromagnetic 
field to subterranean strata; detecting the electromagnetic wave field response; 
seeking, in the wave field response, a component representing a refracted wave; and 
determining the presence and/or nature of any reservoir identified based on the 
presence or absence of a wave component refracted by hydrocarbon layer. 

It further contemplates an apparatus for determining the nature of a subterranean 
reservoir whose approximate geometry and location are known, or for searching for 
a hydrocarbon containing subterranean reservoir, the apparatus comprising: means 
for applying a time varying electromagnetic field to the strata containing the 
reservoir; means for detecting the electromagnetic wave field response; and means 
for seeking, in the wave field response, a component representing a refracted wave, 
thereby enabling the presence and/or nature of a reservoir to be determined. 
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A refracted wave behaves differently, depending on the nature of the stratum in 
which it is propagated. In particular, the propagation losses in hydrocarbon stratum 
are much lower than in a water-bearing stratum while the speed of propagation is 
much higher. Thus, when an oil-bearing reservoir is present, and an EM field is 
5 applied, a strong and rapidly propagated refracted wave can be detected. This may 
therefore indicate the presence of the reservoir or its nature if its presence is already 
known. 

Electromagnetic surveying techniques in themselves are known. However, they are 
10 not widely used in practice. In general, the reservoirs of interest are about 1 km or 
more below the seabed. In order to carry out electromagnetic surveying as a stand 
alone technique in these conditions, with any reasonable degree of resolution, short 
wavelengths are necessary. Unfortunately, such short wavelengths suffer from very 
high attenuation. Long wavelengths do not provide adequate resolution. For these 
15 reasons, seismic techniques are preferred. 



However, while longer wavelengths applied by electromagnetic techniques cannot 
provide sufficient information to provide an accurate indication of the boundaries of 
the various strata, if the geological structure is already known, they can be used to 
20 determine the nature of a particular identified formation, if the possibilities for the 
nature of that formation have significantly differing electromagnetic characteristics. 
The resolution is not particularly important and so longer wavelengths which do not 
suffer from excessive attenuation can be employed. 

* • 

25 The resistivity of seawater is about 0.3 ohm-m and that of the overburden beneath 
the seabed would typically be from 0.3 to 4 ohm-m, for example about 2 ohm-m. 
However, the resistivity of an oil reservoir is likely to be about 20-300 ohm-m. This 
large difference can be exploited using the techniques of the present invention. 
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Typically, the resistivity of a hydrocarbon-bearing formation will be 20 to 300 times 
greater than water-bearing formation. 

Due to the different electromagnetic properties of a gas/oil bearing formation and a 
5 water bearing formation, one can expect a reflection and refraction of the transmitted 
field at the boundary of a gas/oil bearing formation. However, the similarity 
between the properties of the overburden and a reservoir containing water means that 
no reflection or refraction is likely to occur. 

10 Thus, an electric dipole transmitter antenna on or close to the sea floor induces 
electromagnetic (EM) fields and currents in the sea water and in the subsurface 
strata. In the sea water, the EM-fields are strongly attenuated due to the high 
conductivity in the saline environment, whereas the subsurface strata with- less 
conductivity potentially can act as a guide for the EM-fields (less attenuation). If the 

15 frequency is low enough (in the order of 1 Hz), the EM-waves are able to penetrate 
deep into the subsurface, and deeply buried geological layers having higher electrical 
resistivity than the overburden (as e.g. a hydrocarbon filled reservoir) will affect the 
EM-waves. Depending on the angle of incidence and state of polarisation, an EM 
wave incident upon a high resistive layer may excite a ducted (guided) wave mode in 

20 the layer. The ducted mode is propagated laterally along the layer and leaks energy 
back to the overburden and receivers positioned on the sea floor. The term 
"refracted" wave in this specification is intended to refer to this wave mode. 

Both theory and laboratory experiments show that the ducted mode is excited only 
25 for an incident wave with transverse magnetic (TM) polarisation (magnetic field 
perpendicular to the plane of incidence) and at angles of incidence close to the 
Brewster angle and the critical angle (the angle of total reflection). For transverse 
electric (TE) polarisation (electric field perpendicular to the plane of incidence) the 
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ducted mode will not be excited. Since the induced current is proportional to the 
electric field, the current will be parallel to the layer interfaces for TE polarisation 
but, for TM polarisation, there is an appreciable current across the layer interfaces. 

5 A horizontal dipole source on the sea floor will generate both TE and TM waves, but 
by varying the orientation of the receiver antennae, it is possible to vary the 
sensitivity to the two modes of polarisation. It appears that an in-line orientation 
(source and receiver dipoles in-line) is more sensitive to the TM mode of 
polarisation, whereas a parallel orientation (source and receiver dipoles in parallel) is 
10 more sensitive to the TM mode of polarisation. The TM mode is influenced by the 
presence of buried high resistive layers, whereas the TE mode is not. By measuring 
with the two antenna configurations and exploiting the difference between the two 
sets of measurements, it is possible to identify deeply buried high resistivity zones, 
i.e. a hydrocarbon reservoir. 

15 

The present invention has arisen from this realisation. 

According to one aspect of the present invention, there is provided, a method of 
determining the nature of a subterranean reservoir which comprises: deploying an 

20 electric dipole transmitter antenna with its axis generally horizontal; deploying an 
electric dipole receiver antenna in-line with the transmitter; applying an 
electromagnetic (EM) field to the strata containing the reservoir using the 
transmitter; detecting the EM wave field response using the receiver and identifying 
in the response a component representing a refracted wave from the reservoir 

25 according to a first mode; deploying an electric dipole receiver antenna parallel to 
the transmitter; applying an EM field to the strata using the transmitter; detecting the 
EM wave field response using the receiver and identifying in the response a 
component representing a refracted wave from the reservoir according to a second 
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mode; and comparing the first mode refractive wave response with the second mode 
refracted wave response in order to determine the nature of the reservoir. 

According to another aspect of the present invention there is provided, a method of 
5 searching for a hydrocarbon-containing subterranean reservoir which comprises: 
deploying an electric dipole transmitter antenna with its axis generally horizontal; 
deploying an electric dipole receiver antenna in-line with the transmitter; applying an 
EM field to subterranean strata using the transmitter; detecting the EM wave field 
response using the receiver; seeking in the response a component representing a 

10 refracted wave according to a first mode, caused by a high-resistivity zone; 
deploying an electric dipole receiver antenna parallel to the transmitter; applying an 
EM field to the strata using the transmitter; detecting the EM wave field response 
using the receiver; seeking in the response a component representing a refracted 
wave according to a second mode; and comparing the first mode refractive wave 

15 response with the second mode refractive wave response in order to determine the 
presence and/or nature of any high-resistivity zone. 

The first mode may be considered to be a TM mode, and the second mode a TE 
mode. 

20 

Thus, according to the invention, measurements are taken with the transmitter and 
receiver both in-line and parallel, and the two sets of measurements are compared. 
A characteristic difference in values indicates a highly resistive layer located beneath 
highly conductive strata. High resistivity indicates the presence of hydrocarbons and 
25 so the difference in values is a direct hydrocarbon indicator. 

This technique can be used in conjunction with conventional seismic techniques to 
identify hydrocarbon reservoirs. 
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Preferably, the transmitter and/or receiver comprises an array of dipole antennae. 

The technique is applicable in exploring land-based subterranean reservoirs but is 
5 especially applicable to submarine, in particular sub-sea, subterranean reservoirs. 
Preferably the field is applied using one or more transmitters located on the earth's 
surface, and the detection is carried out by one or more receivers located on the 
earth's surface. In a preferred application, the transmitter(s) and/or receivers are 
located on or close to the seabed or the bed of some other area of water. 

10 

In a preferred arrangement, the transmitter and receiver antennae are located on a 
common cable towed behind a vessel. This will result in a fixed offset or a series of 
fixed offsets where several receivers are employed. Preferably, the transmitter 
transmits both modes and may therefore comprise two dipoles arranged mutually at 
1 5 right angles . Preferably each receiver comprises two dipoles mutually at right 

angles. Preferably one transmitter dipole and one receiver dipole are arranged at 
right angles to the direction of the cable. Alternatively the transmitter and/or 
receivers may each comprise a single dipole antenna arranged obliquely, eg. at 45° 
to the direction of the cable. With this arrangement the transmitted field is resolved. 

20 

Using this technique, it is possible to achieve comparable results from the two modes 
as the same signal and offset are used. It will not matter greatly if the transmitter 
drifts in frequency or amplitude. Furthermore, reservoirs can be detected in real 
time. Thus, if the results show a difference in the two modes, this will strongly 
25 indicate the presence of an H/C bearing reservoir and so a more detailed study can 
be made at once. 

Such a system would generally use a single transmission source and several 
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receivers, typically more than ten. The different offsets would be suitable for 
detecting reservoirs at different depths. 

The receivers can be deployed on a single cable or on a series of parallel cables. 
5 There may also be several transmitters. 

In practice, the vessel would normally stop and the cable allowed to sink prior to 
transmission. There would be a transmission at several different frequencies before 
moving to another location. The technique is particularly suitable for edge detection, 
10 and it is a simple matter to select a suitable resolution. However, if the surveying is 
being carried out in an undetermined area, the resistivity in the top layers should be 
mapped, for example with MT methods or by inversion after a reflection study. 

If the offset between the transmitter and receiver is significantly greater than three 
15 times the depth of the reservoir from the seabed (i.e. the thickness of the 

overburden), it will be appreciated that the attenuation of the refracted wave will 
often be less than that of direct wave and the reflected wave. The reason for this is 
the fact that the path of the refracted wave will be effectively distance from the 
transmitter down to the reservoir i.e. the thickness of the overburden, plus the offset 
20 along the reservoir, plus the distance from the reservoir up to the receivers i.e. once 
again the thickness of the overburden. 

The polarization of the source transmission will determine how much energy is 
transmitted into the oil-bearing layer in the direction of the receiver. A dipole 
25 antenna is therefore the selected transmitter. In general, it is preferable to adopt a 
dipole with a large effective length. The transmitter dipole may therefore be 100 to 
1000 meters in length and may be towed in two orthogonal directions. The receiver 
dipole optimum length is determined by the thickness of the overburden. 
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The transmitted field may be pulsed, however, a coherent continuous wave with 
stepped frequencies is preferred. It may be transmitted for a significant period of 
time, during which the transmitter should preferably be stationary (although it could 
5 be moving slowly), and the transmission stable. Thus, the field may be transmitted 
for a period of time from 3 seconds to 60 minutes, preferably from 3 to 30 minutes, 
for example about 20 minutes. The receivers may also be arranged to detect a direct 
wave and a wave refracted from the reservoir, and the analysis may include 
extracting phase and amplitude data of the refracted wave from corresponding data 
1 0 from the direct wave. 

Preferably, the wavelength of the transmission should be in the range 
0.1s<X<5s; 

15 

where X is the wavelength of the transmission through the overburden and s is the 
distance from the seabed to the reservoir. More preferably X is from about 0.5s to 
2s. The transmission frequency may be from 0.01 Hz to 1 kHz, preferably from 1 
to 20 Hz, for example 5 Hz. 

20 

Preferably, the distance between the transmitter and a receiver should be in the range 
0.5?i<L<10V, 

» ■ 

25 where X is the wavelength of the transmission through the overburden and L is the 
distance between the transmitter and the first receiver. 

It will be appreciated that the present invention may be used to determine the position, 
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the extent, the nature and the volume of a particular stratum, and may also be used to 
detect changes in these parameters over a period of time. 

The present invention also extends to a method of surveying subterranean measures 
5 which comprises; performing a seismic survey to determine the geological structure of 
a region; and where that survey reveals the presence of a subterranean reservoir, 
subsequently performing a method as described above. 

The invention may be carried into practice in various ways and will now be illustrated 
10 in the following embodiments and reduced scale investigations and simulations. In the 
accompanying drawings, 

Figure 1 is a vertical cross-section through a testing tank; 
1 5 Figure 2 is a plan view of the tank of Figure 1 ; 

* 

Figure 3 is a plan view of the antennae used in the tank of Figure 1 ; 
Figure 4 is a side view of the antenna in Figure 3 ; 

20 

Figures 5 and 6 are respectively a schematic plan view and side view of the testing 
tank set up for measurement; 

Figures 7 is a graph showing calculated and measured values for the transmitted 
25 electric field for a given frequency in the model experiment; 

Figure 8 is a graph showing calculated values for the electric field in a realistic earth 
model; 

30 Figure 9 is a schematic side view of a cable layout towed by a vessel; 
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Figure 10 is apian view corresponding to Figure 9; and 

Figures 11 and 12 are views similar to Figure 10 showing two alternative 
5 arrangements. 

The tank 1 1 shown in Figure 1 and 2 comprises a concrete enclosure 9m long, 6m 
wide and 8m in depth. The tank 1 1 is filled with sea water 12. A diaphragm 13 filled 
with fresh water 14 is located in the tank. The diaphragm 13 is 7.5m long, 4.25m wide 
1 0 and 0.25m thick and can be located at any desired height in a horizontal orientation 
within the tank 11. 

The conductivity of the sea water 12 was measured to be 5.3 S/m at 14°C and the 
conductivity of the fresh water was measured to be 0.01 3 S/m. The ratio of the two 
1 5 conductivities is therefore very close to 400. 

The critical frequency f c of a conducting medium, i.e. the frequency at which the 
displacement current is equal to the conduction current, is given by 

20 fi = — — »18— GHz 

IftSrSa Sr 

where s r is the relative dielectric constant of the medium, and a the conductivity in 
S/m. For water, s r - 80 at the frequencies and temperatures of interest. For the two 
25 conductivity values a - 5.2 S/m and a = 0.013 S/m, f c = 1.2 GHz and 3 Mhz, 

respectively. Since, in the experiments, the highest frequency is 0.83 MHz, it is a fair 
approximation to neglect the displacement current, even for the fresh water. 

For a non magnetic, conductive medium, the propagation constant y is given by 
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If & ^ 

O) \ SrSo + — i 

\K J<*>) 



y = Ja>y[sju — ja> I srSo + — — juo 



^ ]G)f.io<y — lir^^-- 



fir 



The wavelength \ defined as the distance in which the phase changes 2tt 3 is given by 



X 2 * 



Re{r} 



X m,/in MHz and a in S/m. The skin depth, the distance in which the amplitude 
diminishes by 1/e, is related to the wavelength by 



71 



2fjUo<T 



For the extremes of the frequency range, for the sea water with a = 5.2 S/m, 



Frequency 30 kHz 830 kHz 

Skin depth 1 .27 m 0.24 m 

Wavelength 8.01 m 1 .52 m 



and for the fresh water with o = 0.013 S/m, 

Frequency 30 kHz 830 kHz 

Skin depth 25.4 m 4.8 m 

Wavelength 160.2 m 30.4 m 
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Referring now to Figures 3 and 4, two identical electrical dipole antennae, as shown 
were used for the transmitter and receiver. 

Each antenna 15 comprises two square brass plates 16, 15cm square, mounted on an 
5 epoxy substrate 1 7. Each plate 1 6 is connected to a co-axial cable 1 8, which passes 
through an epoxy tube 19 mounted at right angles to the plate 16, to a balun which 
transforms the impedance of the antenna 15 from about 2£2 in sea water to about 5011. 

The measurement set-up is shown in Figure 5 and 6. An automatic network analyser 
1 0 (ANA) measures the transmission between the antennae 15 as a function of distance 
(offset) and frequency. The arrangement shown in Figure 5 shows the antennae 15 in 
the parallel orientation. The in-line orientation is achieved by rotating both antennae 
through 90° in the horizontal plane. 

1 5 The results of the measurements are shown in Figure 7 together with the 

corresponding theoretical results. The measurements agree well with the theoretical 
results and the figure contains two sets of curves, one with parallel antennae and one 
with the antennae in line. The theoretical results are computed for infinitesimal dipole 
antennae. The orientations of the antennae and the frequency are shown on the 

20 Figures. 

The parameters of the experiment are scaled relative to possible practical situations. 
To give an idea of orders of magnitude; if the frequency is scaled down by a factor of 
40,000 and the conductivity by a factor of 10, the dimensions will be scaled up by a 
25 factor of 632, and the experimental setup would correspond to a low conductivity layer 
of thickness 150 m and conductivity 0.0013 S/m below an overburden of thickness 
300 m and conductivity 0.52 S/m. The corresponding frequency range would be from 
0.75 Hz to 20 Hz, and the length of the antenna nearly 300 m. 

a 

30 The method with the TE and TM-mode have been tested by computer simulations on a 
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simple horizontally layered earth model with electrical parameter values for typical 
deep water subsurface sediments. The model has an infinite insulating air layer, a 
1150 meter water layer of 0.3 1 25Qm, 950 meter overburden of 1 Clm, a 1 50 meter 
reservoir zone of 50£lm and an infinite underbidden of lQm. Figure 8 illustrates the 
5 amplitude response | E | (electric field) as a function of receiver offset, caused by a 1 
Hz signal. Responses from both the TM-mode (solid with x f s) and TE-mode (dashed 
with +'s) are shown. The amplitudes for the TM-mode are approximately 10 times 
larger at an offset of 5 km. As a reference, the response from a homogenous half- 
space of 1 Qm is shown for both configurations (corresponding to a response from a 
10 water filler reservoir or outside the reservoir area). The TE-mode has the largest 

deviation from its half-space, ie. this mode is more sensitive to a hydrocarbon layer. 

Figures 9 and 1 0 show a vessel 3 1 towing a cable (or streamer) 32 just above the 
seabed 33. The cable 32 carries a transmitter dipole antenna 34 and several receiver 
1 5 dipoles 35, only four of which are shown. The depth of water might be of the order of 
1 000m, the offset between the transmitter 34 and the nearest receiver 35 might be 
about 2000m and the receivers might be about 100m apart. The transmitter 34 is 
controlled from the vessel 3 1 via the cable 32 and the responses detected by the 
receivers 35 are relayed back to the vessel 3 1 in real time, again via the cable 32. 

20 

Figure 10 shows an arrangement in which the vessel 31 tows three cables 41, 42, 43, 
each carrying a series of receivers 45, 46, 47. The spacing of the three cables 41, 42, 
43 is achieved by means of a spar 44. 

25 In the arrangement shown in Figure 1 1 , the transmitter 48 is in the form of two dipole 
antennae, one parallel to the cable 42 and one at right angles. 

The arrangement shown in Figure 12 is similar to Figure 1 1, but in this case, the 
transmitter 5 1 is a single dipole antenna arranged at 45° to the cable 42. 
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Claims 

1. A method of determining the nature of a subterranean reservoir which 
comprises: deploying an electric dipole transmitter antenna with its axis generally 

5 horizontal; deploying an electric dipole receiver antenna in-line with the transmitter; 
applying an electromagnetic (EM) field to the strata containing the reservoir using 
the transmitter; detecting the EM wave field response using the receiver and 
identifying in the response a component representing a refracted wave from the 
reservoir according to a first mode; deploying an electric dipole receiver antenna 

10 parallel to the transmitter; applying an EM field to the strata using the transmitter; 
detecting the EM wave field response using the receiver and identifying in the 
response a component representing a refracted wave from the reservoir according to 
a second mode; and comparing the first mode refractive wave response with the 
second mode refracted wave response in order to determine the nature of the 

15 reservoir. 

2. A method of searching for a hydrocarbon-containing subterranean reservoir 
which comprises: deploying an electric dipole transmitter antenna with its axis 
generally horizontal; deploying an electric dipole receiver antenna in-line with the 

20 transmitter; applying an EM field to subterranean strata using the transmitter; 
detecting the EM wave field response using the receiver; seeking in the response a 
component representing a refracted wave according to a first mode, caused by a 
high-resistivity zone; deploying an electric dipole receiver antenna parallel to the 
transmitter; applying an EM field to the strata using the transmitter; detecting the 

25 EM wave field response using the receiver; seeking in the response a component 
representing a_ refracted wave according to a second mode; and comparing the first 
mode refractive wave response with the second mode refractive wave response in 
order to determine the presence and/or nature of any high-resistivity zone. 
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3. A method as claimed in Claim 1 or Claim 2, characterised in that the first 
mode is a TM mode of polarisation and/or the second mode is a TE mode of 
polarisation. 

5 

4. A method as claimed in any preceding Claim, characterised in that the 
transmitter and/or receiver comprises an array of dipole antennae. 

5. A method as claimed in any preceding Claim, characterised in that the 
10 transmitter and/or receiver is located on or close to the seabed or the bed of some 

other area of water. 

6. A method as claimed in any preceding Claim, characterised in that the 
transmitter and receivers are located on a common cable arranged to be towed 

15 behind a vessel. 

7. A method as claimed in any preceding Claim, characterised in that the 
transmitter comprises two dipole antennae arranged mutually at right angles. 

20 8. A method as claimed in any preceding Claim, characterised in that each 
receiver comprises two dipole antennae arranged mutually at right angles. 

9. A method as claimed in Claim 6 characterised in that the transmitter and/or 
receiver each comprise a single dipole antenna arranged obliquely to the direction of 

25 the cable. 

10. A method as claimed in any preceding Claim, characterised in that the 
frequency of the EM field is continuously varied over the transmission period. 



■ 
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11. A method as claimed in any preceding Claim, characterised in that the field is 
transmitted for a period of time for 3 seconds to 60 minutes. 

5 12. A method as claimed in Claim 11, characterised in that the transmission time 
is from 3 to 30 minutes. 

13. A method as claimed in any preceding Claim, characterised in that the 
wavelength of the transmission is given by the formula 

0.1s < X <10s; 

■i 

wherein X is the wavelength of the transmission through the overburden and s is the 
distance from the seabed to the reservoir. 

14. A method as claimed in any preceding Claim, characterised in that distance 
between the transmitter and a receiver is given by the formula 

0.5 < A, < 10s; 
20 

where X is the wavelength of the transmission through the overburden and 1 is the 
distance between the transmitter and the receiver. 

15. A method as claimed in any of Claims 10 to 14, characterised in that the 
25 transmission frequency is from 0.01 Hz to 1 kHz. 

16. A method as claimed in Claim 15, characterised in that the transmission 
frequency is from 1 to 20 Hz. 
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17. A method as claimed in any preceding Claim, characterised in that it includes 
suppressing the direct wave and/or any other known wave contribution that may 
disturb the measurements, thereby reducing the required dynamic range of the 
receiver and increasing the resolution of the refracted wave. 

18. A method of surveying subterranean measures which comprises: performing 
a seismic survey to determine the geological structure of a region and where that 
survey reveals the presence of a subterranean reservoir, subsequently performing a 
method as claimed in any preceding Claim. 
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